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^ ' Abstract. - We investigate the efi'ect of surface waves generated by the Faraday instabihty on a 

, shear-thickening surfactant solution under vertical vibrations. We show that a prolonged oscillation 

of the surface above the instability onset leads to an increase of the fluid viscosity. This phenomenon 
is evidenced by comparing the time needed for the instability to develop in the fluid at rest and 
, once the fluid has been shaken above onset: pre-shaking may delay the instability by two orders 

(~] ' of magnitude. A simple model based on a time-dependent viscosity is proposed which accounts 

Q ' quantitatively for the experimental observations. 

o ■ 



> 
o 

Introduction. — In the last few decades, the coupling between microstructure and 
I flow in complex fluids has been the subject of intense research effort [1]. In most cases the 

. ■ flow is a simple shear flow induced by a rheometer. Yet, another way to induce a controlled 

I flow is by the means of hydrodynamic instabilities. In this paper the effect of the well- 



known Faraday instability, by which the initially flat surface of a fluid layer submitted to 
I vertical vibrations gives way to surface waves with characteristic wave number kc above 

some critical acceleration ac [2], is tested on a shear-thickening dilute solution of wormlike 
micelles [3] . Recently the Faraday instability in complex fluids has been the subject of various 
■ studies [4-11]. In most of these previous works [4-8], the authors compared the behaviour 

H I of a complex fluid to that of a Newtonian one close to the instability threshold. Differences 

were explained by invoking the linear viscoelastic properties of the samples. Far above the 
instability threshold, a stronger coupling between the flow induced by the instability and 
the complex fluid microstructure is expected. For instance, in clay suspensions, hysteresis 
and high amplitude fingerlike structures were shown to be linked to shear-thinning [9]. In 
shear-thickening colloidal suspensions, external perturbations may grow and get stabilized 
as "persistent holes" across the fluid surface [10]. In a dilute polymer solution, the flow- 
microstructure coupling can delay the transition to disordered states [11]. The present work 
shows that, far above the instability threshold, the surface waves become large enough to 
induce shear-thickening in dilute wormlike micelles. In the following we first present the 
complex fiuid under study and the experimental setup. Then measurements of the set-up 
time of the Faraday instability, defined as the time needed for the surface waves to become 
large enough to be detected by our experimental setup, are presented. We show that the 
set-up time may increase by two orders of magnitude when the fluid is pre-shaken above the 
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Fig. 1: (a) Temporal evolution of the viscosity of a 250 ppm CTAB-NaSal solution for 7 = 5 s ^ 
at T = 20° C. (b) Viscosity 7; as a function of the shear rate 7. See text for experimental details. 

instability threshold. Finally a model is proposed to account for such observations, based 
on a time-dependence of the viscosity induced by the surface waves. 

System under study and rheological measurements. Our working fluid is a di- 
lute solution of wormlike micelles made of CTAB-NaSal in water at a very low concentration 
ccTAB = 250 ppm by weight and the molar ratio of salt to surfactant is 1:1. This wormlike 
micelle solution in the dilute regime is known to display a shear-thickening behaviour at 
typical shear rates 7 = 1-10 s~^ [3]. Such shear-thickening has also been observed in a 
wide range of other wormlike micelle solutions and appears as a common feature of dilute 
systems [12-14]. In order to characterise our low- viscosity micellar system, rheological mea- 
surements are performed in a double Couette geometry (rotor radii 28 and 34 mm, stator 
radii 29.5 and 32 mm) thermostated at T = 20°C using a strain-controlled Rheometrics 
RFSII rheometer. Figure 1(a) shows the dynamic viscosity 7y as a function of time when a 
shear rate of 5 s~^ is imposed at time t = 0. After an induction time of about 200 s, rj^t) 
increases by a factor of 10 within a characteristic time of about 500 s. This temporal evolu- 
tion of the viscosity is comparable to that reported in the same system by Liu and Pine [3] . 
Next a constitutive curve 77 vs 7 is obtained by submitting our fluid to various shear rates 
7 for 1000 s and by averaging r]{t) over the last 100 s. Such a procedure is repeated three 
times. The viscosity shown in fig. 1(b) is the mean of these three measurements and the 
error bars are their standard deviation. It is clearly seen that between 7 = 0.5 and 5 s""'^, the 
viscosity increases with increasing shear rate up to = 16 mPa.s. Such a shear-thickening 
phenomenon was attributed to the formation of a shear-induced gel-like structure [15,16], 
which was confirmed by light scattering experiments [12]. At greater shear rates, the fluid 
becomes shear-thinning as expected for a gel that gets progressively broken down by shear. 
The curve 77(7) of fig. 1(b) is only qualitatively similar to that of ref. [3]. In our opinion, the 
large error bars shown by our data in the shear-thickening regime illustrate the variability 
in both the final value of the viscosity and in the induction and growth times of ri(t), as well 
as the sensitivity of viscosity measurements to filling conditions and contamination in such 
dilute samples. 

Experimental setup for the Faraday instability. Our experimental setup has 
been described in ref. [8]. A cylindrical cell of depth 10 mm and diameter 60 mm filled with 
our fluid and thermostated at T = 20±0.1°C is submitted to vertical vibrations of frequency 
/ by an electromagnetic shaker (Ling Dynamic Systems V406). The fluid layer is lit from 
above so that the surface is seen bright only when it makes an angle greater than 1° with 
the horizontal plane. Close to the instability onset, small amplitude surface patterns can 
be approximated by a superposition of sinusoidal waves of wave number kc and amplitude 
^. With our setup, surface waves are detected only when ^ > = tan(7r/180)/fcc. In our 
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experiments, fee — 1000 m^^ so that ~ 20 /im. 

In a first step, measurements of the critical acceleration and wave number (not 
shown) were performed. As reported earlier by Raynal et al. [4], the viscoclastic behaviour 
of the solution leads to Oc values lower than those of a Newtonian fluid with the same 
zero-shear viscosity while kc follows the same dispersion relation. Besides, in our system, 
the response of the surface is always subharmonic [i.e. at f/2) over the frequency range 
/ = 20-200 Hz and Uc measurements do not reveal any significant hysteresis. Two elements 
may explain the absence of any signature of shear-thickening in these measurements. First 
the time needed to measure Uc and kc at a given frequency (~ 3 min) is smaller than the 
characteristic time for shear-thickening (~ 10 min). Second although the equivalent shear 
rate, estimated as 7 = Tr/^^^c close to instability onset [6], falls into the shear-thickening 
regime (one finds 7 ~ 6 s~^ for / = 100 Hz), the wave amplitude may still be too small to 
induce any modification of the fluid microstructure because the flow involved in the Faraday 
experiment is both oscillating and elongational and thus strongly differs from the continuous 
shear used in rheological experiments. Therefore, in order to induce shear-thickening through 
the Faraday instability, the experiments detailed below were performed on longer time scales 
and at accelerations far above the instability threshold where the wave amplitude is larger. 

Experimental results. — To evidence a signature of shear-thickening, we measure 
the set-up time ta of the instability in our fluid once it has been shaken at an acceleration 
Uy above onset for a duration ty. The exact experimental procedure is sketched in fig. 2(a). 
The fluid layer is first vibrated at frequency / and ay = 2ac for a duration ty. Such a 
"pre-shaking" is followed by a short period when the instability is stopped (a = ac/2 for 3 s 
then a = Q.9ac for 1 s) and the surface wave completely disappears. Finally the acceleration 
is set to a = (1-1- ejUc and the set-up time, i.e. the time needed for the surface wave to 
be detected again by our experimental set-up, is measured. This method is preferred to a 
direct measurement of Oc because the time needed to precisely measure Qc may vary from 
1 to 5 min and therefore is not suitable for a time-dependent microstructure. In all cases, 
two successive measurements of ta are separated by at least ten minutes, in order to make 
sure that our fluid is not perturbed by previous excitation. 

In a first experiment, ty is varied from 1 to 1800 s for / = 100 Hz and e = 0.05. Fig- 
ure 2(b) shows that ta strongly increases with increasing ty, up to two orders of magnitude. 
More precisely, ta remains roughly constant up to ty ~ 200 s which corresponds to the 
induction time mentioned above in fig. 1(a). ta then saturates to about 600 s which again is 
consistent with the characteristic growth time of the shear-induced structure. We interpret 
this remarkable effect on the set-up time as a consequence of shear-thickening induced by 
pre-shaking above onset. 

In a second experiment, ta is measured after pre-shaking at fixed Uy = 2ac and ty = 600 s 
for various final accelerations a above Oc. The results are shown in fig. 3 and compared to the 
set-up time measured in the fluid initially at rest {i.e. without pre-shaking) for three different 
frequencies. Let us first discuss the results obtained in the fluid at rest. As expected, one 
finds ta oc 1/e, where e = a/ac — 1 is the dimensionless distance to the threshold. Indeed, 
in general, as long as the amplitude ^ of the surface wave is small enough, £^{t) follows the 
equation [17]: 

(1) 

where e/xg is the growth rate, so that for a perturbation of initial amplitude ^Oi the set-up 
time ta deflned by ^(ta) = S.d is given by: 

ta^To^^, (2) 
a — Oc 

where Tq = Tg hj{^d/^o) is the set-up time for e = 1 (i.e. a = 2ac). As shown in fig. 3, eq. (2) 
perfectly matches the experimental data in the fluid at rest and allows us to determine 
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Fig. 2: (a) Experimental procedure, (b) Set-up time ta versus the pre-shaking duration tv at 
a„ = 2ac and / = 100 Hz. The continuous line is the result of the best fit by eqs. (8) and (9) with 
T = 461 s and 6ao = 0.149. 



precisely Cc hereafter noted uq for the sake of clarity (to within less than 1%) and Tq (to 
about 2%). The best fit parameters are given in table 1. 

Let us now turn to the set-up time in the pre-shaken fluid. As seen in fig. 3, for a given 
e, ta is much larger (sometimes a hundred times larger) for the pre-shaken fluid than for 
the fluid at rest. In the following, we argue that this behaviour is linked to an increase of 
the fluid viscosity induced by the surface waves when the fluid is shaken at a„ = 2ao. After 
the instability is stopped and the acceleration is set to its final value, the viscosity tends 
to relax towards that of the fluid at rest. In this case, the instability threshold and the 
characteristic time Tg become time-dependent. The set-up time thus results from a trade-off 
between the slow relaxation of the viscosity and the increase of the critical acceleration and 
growth rate induced by shear-thickening. In order to get a more quantitative understanding 
of this phenomenon, a simple model is proposed thereafter. 

Model. — Our main assumption is that, once the surface waves have disappeared and 
the acceleration is set to a at time t = 0, the viscosity decreases exponentially with a 
characteristic time r: for t > 0, i^{t) = 7;o(l -I- (5exp(— t/r)), where t/o is the viscosity of 
the fluid at rest and 5 accounts for shear-thickening induced during pre-shaking so that 
i]{t = 0) = 770(1 + 5). Moreover, it is known that, for a Newtonian fluid, Oc oc 77 and 
Tg cx rj^^ [17]. If we assume our dilute micellar solution to be a Newtonian fluid whose only 
time-dependent parameter is the viscosity, we have: 

ac{t) = ao— =ao(l + ^e-*/^), (3) 
'70 

where e/rgo is the growth rate for the fluid at rest. 

For a given acceleration a above oq, two cases arc possible: a < adO) or a > ac(0). If 
a < ac(0), the surface remains flat as long as a < ac{t). Let ti be the time needed for the 
viscosity to relax so that a = adti). Equation (3) yields ti = rln((5/e). By changing the 
time origin t ^ t — ti, we are left with solving the case a = 0^(0) i.e. e = 6. For small 
amplitudes and time-dependent Tg and e, eq. (1) becomes 



dt \ac{t) 



Inserting eqs. (3) and (4) and integrating eq. (5) over time, one flnds To/e = ^2 — '''(I — 
exp(— t2/T)), where t2 is the set-up time in the case a ~ ac(0) and Tg = Tgo ln(^d/Co)- If one 
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a (m.s ) 
7.37 8.02 



a (m.s ) 
10.98 11.96 



a (m.s ) 
19.46 21.2 




Fig. 3: Set-up time ta versus a (top scale) and e (bottom scale) for the fluid at rest (•) and for the 
pre-shaken fluid (o). The continuous line is the best fit by eqs. (6) and (7). The vibration frequency 
is (a) / = 60 Hz, (b) / — 85 Hz, and (c) / = 135 Hz. The flt parameters are given in table 1. 



further assumes that t2/T <^ 1, an expansion to second order in t2/T leads to t2 ~ jt. 
Finally, when a < ac(0), the set-up time ta is the sum of the two terms ti and t2'. 



ta = tl + t-I = T hi 



2ToT 



(6) 



The first term corresponds to the decrease of ac{t) down to the imposed acceleration, while 
the second term derives from the growth of the surface wave amplitude up to ^d- 

Turning now to the case a > ac(0), one has e{t) = (e — ^e^*/"^)/ {6e^*^'^ + 1). Integrating 
eq. (5) under the assumption ia/r <C 1 leads to 



tn 



1 



1 



2ToT 



(7) 



It is easily checked that for e = S, eqs. (6) and (7) give the same value for ta and that for 
S, the same scaling as for the fluid at rest, namely ta = 7o/e, is recovered from eq. (7). 
Although our expression for ta involves four parameters, uq, Tq, S, and r, there are 
actually only two free parameters S and r. Indeed, as already mentioned, uq and Tq are 
determined experimentally by fitting ta(e) for the fiuid initially at rest using eq. (2). Then 
6 and r are found by fitting ta{e) for the pre-shaken fluid with eqs. (6) and (7), where ao 
and To are fixed to the values found previously. The results of these fits are shown in fig. 3 
and the various fit parameters are presented in table 1. The fits are in good quantitative 
agreement with the experimental results. As expected, increases as / increases and r 
remains roughly constant r ~ 520 s ±10%. A posteriori, one can check that t2 t a.s soon 
as e > 0.002. 



/ (Hz) 


65 


85 


135 


ao (m.s^^) 


7.3 


10.9 


19.3 


To (s) 


0.49 


0.35 


0.28 


r (s) 


501 


565 


502 


S 


0.136 


0.105 


0.076 



Table 1: Best fit parameters used in eqs. (6) and (7) to model the ta vs e data for / = 65, 85, and 
135 Hz. The resulting fits are shown in flg. 3. 
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Going back to the first experiment reported in fig. 2(b), the same arguments can be used 
to understand the increase of ta with ty at a fixed e. Here wc assume that during pre-shaking 
at Uv = 2ao, the viscosity increases with ty as r]{tv) = ?7o(l + i5oo(l — exp(ti,/r))), where 
77o(l + <5oo) is the asymptotic viscosity and r the characteristic time of shear-thickening. Note 
that, although the exponential decay proposed above is probably a reasonable model for the 
relaxation of the viscosity towards its value at rest, an exponential growth of 7/(t) is clearly 
inconsistent with the rheological measurements shown in fig. 1(a) which rather involve an 
induction time of about 200 s. In the present approach, we chose to keep an exponential 
form for 77 (t) so that the equations remain tractable analytically. Following the path used to 
find eqs. (6) and (7), we end up with the following expressions: if Soo{l — e^*"/'^) > e = 0.05, 

>.^^2.r„(l + l)+.ln( ^-"7"-^^' ), (8) 

and if daoil - e-*"/^) < e = 0.05, 

^^{'- 5^(1 - e-^V.)) + f ^ (1 - ^^(1 _V^^V.)) + sZiT^^y 

As shown in fig. 2(b), the fit obtained with a characteristic time t = 461 s and Soo = 0.149 
closely matches the experimental data. Such an accuracy of our crude exponential model 
for ri{t) points to a weak dependence of the set-up time on the details of rj{t). 

Discussion. — Let us now discuss the present results and model. The various relaxation 
times T found by fitting our simple model to the experimental data are strikingly similar 
to the characteristic times revealed in rheological experiments. We conclude that the same 
shear-thickening phenomenon is at play under pure shear flow and in the surface wave 
pattern induced by the Faraday instability. As in previous work under shear [12], we tried to 
evidence turbid regions characteristic of a gel-like structure using light scattering. However 
such experiments did not reveal any significant increase of the turbidity neither locally nor 
globally. This is probably because the effect is too small or too localized compared to the 
one at constant shear rate. Indeed, as indicated by our model, the viscosity increase is 
much smaller under vibrations {5 ^ 10%) than under continuous shear (~ 1000%). We 
believe that the main reason for this discrepancy is that, in the Faraday instability, the 
shear rate at the surface oscillates in time and the period of the oscillation is much smaller 
than T (1// = 0.01 ^ r ~ 500 s). Therefore, one should rather compare the shear strains 
T7 ~ 2000 under continuous shear to iir^kc ~ 2 in the Faraday experiment, where 7 ~ 5 s~^ 
is the typical shear rate for shear-thickening and ^ ~ 300 /im is the amplitude of the surface 
waves at e = 1 . The fact that T7 ^ could explain why the increase of viscosity is much 

smaller in the Faraday experiment but still noticeable through set-up time measurements. 
Finally, note that the comparison with pure shear flow (which involoves both rotation and 
elongation) may not be that relevant since the flow field induced by the Faraday instability is 
purely elongational besides a viscous sublayer [11]. In the absence of any published data on 
the present micellar system under elongational flow, the issue of the exact interplay between 
Faraday waves and shear-thickening is left as an open question. 

Conclusion. — Set-up time measurements of the Faraday instability revealed a shear- 
thickening phenomenon induced by surface waves. A simple phenomenological model was 
proposed to account for the experimental observations. Even if no gel-like structure was 
directly evidenced, these experiments show that, far above the instability threshold, the 
surface deformation is large enough to induce a change in the microstructure of this partic- 
ular complex fluid. Further experiments will deal with fluids whose shorter relaxation times 
may couple more strongly to the Faraday instability. An interesting question would then be 
whether the structuration of the flow into patterns with some characteristic wave number 
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can be used to imprint some macroscopic spatial organization of the fluid microstructure 
with the same wave number. 
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